Time-dependent control of ultracold atoms and Bose-Einstein condensates in 

magnetic traps 

N. V. Vitanov 1 and K.-A. Suominen 1,2 
(1) Helsinki Institute of Physics, PL 9, FIN-00014 Helsingin yliopisto, Finland 
(2) Theoretical Physics Division, Department of Physics, University of Helsinki, PL 9, FIN-00014 Helsingin yliopisto, Finland 



o 



O 



> ■ 
in ■ 

<n : 
o . 
r» ■ 
o ■ 

: 

On . 

o : 

>>: 
43 : 



- 1—^ 

X 



With radiofrequency fields one can control ultracold atoms 
in magnetic traps. These fields couple the atomic spin states, 
and are used in evaporative cooling which can lead to Bose- 
Einstein condensation in the atom cloud. Also, they can be 
used in controlled release of the condensate from the trap, 
thus providing output couplers for atom lasers. In this pa- 
per we show that the time-dependent multistate models for 
these processes have exact solutions which are polynomials of 
the solutions of the corresponding two-state models. This 
allows simple, in many cases analytic descriptions for the 
time-dependent control of the magnetically trapped ultracold 
atoms. 
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Neutral atoms possessing hyperfme structure can be 
trapped in spatially inhomogeneous magnetic fields, if 
they are in the appropriate spin state. As the magnetic 
field imposes spin-dependent Zeeman shifts on the atomic 
energy levels, a spatially changing magnetic field maps 
into an external potential felt by the atom (see Fig. [I]). 
The magnetic traps are, however, very shallow, so only 
ultracold atoms can be trapped. For alkali atoms the 
proper temperatures have been obtained via precooling 
with laser light. Once the atoms are trapped, one can 
decrease the trap depth, which allows the hot atoms to 
escape, and those left behind thermalize via collisions 
into a lower temperature — this is called evaporative cool- 
ing Q. By combining magnetic traps with evaporative 
cooling one can now reach the densities and temperatures 
where Bose-Einstein condensation takes place [0. Then 
the atoms form a coherent superposition, which can be 
released from the trap j^j. As the escaping atoms main- 
tain their coherence Q] , the experiment is a prototype for 
an atom laser, i.e., production of coherent, propagating 
packets of matter waves. 

Evaporative cooling requires effective and precise con- 
trol of the trap depth. This can be achieved by cou- 
pling the spin states (labeled with the magnetic quantum 
number m) with a radiofrequency (rf) field Jl],||||. The 
coupling introduces transitions between the neighbouring 
spin states (Am = ±1). As demonstrated in Fig. |l|, the 
frequency cj r f of the field controls the depth of the trap. If 
we want to estimate the efficiency of evaporative cooling, 
we can transform the rf resonances into curve crossings, 
and describe the dynamics of the atoms at the edge of the 
trap with an appropriate time-dependent curve crossing 
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FIG. 1. Magnetic trapping of neutral atoms. We show 
the atomic potentials for the / = 1 trapping system, and 
how the radiofrequency field (frequency ui T f) can be used 
to control the trap depth fcsTd. The resonant transition 
m = —1 — » m = — » m = +1 moves the atoms with ki- 
netic energy larger than ksTd from the trapping state to a 
strongly nontrapping state. 

model jj], see Fig. [| Typically one compares the trap os- 
cillation period times the spin-change probability at the 
resonance point to the other time scales of the trapping 
and cooling process, including collisional loss rates 0. 

However, as one usually operates in the region where 
Zeeman shifts are linear, the rf field couples sequentially 
all spin states, instead of just selecting a certain pair of 
states. In case of only two spin states we can use the stan- 
dard Landau-Zener model in estimates of the efficiency of 
evaporative cooling. In practice, however, one has 2/ + 1 
states, where / is the hyperfme quantum number of the 
atomic state used for trapping. So far condensates have 
been realised for f — 1 and / = 2, but experiments 
e.g. with cesium involve states with / = 3 and / = 4, so 
in order to achieve efficient evaporative cooling one needs 
several sequential rf- induced transitions. Clearly the use 
of the two-state Landau-Zener model can be questioned 
in these multistate cases. 

Once the Bose-Einstein condensation has been 
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achieved one can release the condensate just by switch- 
ing the magnetic field off. This technique, however, does 
not allow much control over the release. Moreover, it al- 
ways involves all the atoms. With rf fields one can trans- 
form parts of the condensate into untrapped states, in 
which they are typically accelerated away — this has been 
demonstrated experimentally || , as well as the fact that 
the released atoms are in a coherent superposition |Q, 
thus justifying the term "atom laser" . This output cou- 
pling process can be achieved either by using rf pulses 
which are resonant at the trap center, or by using chirped 
rf fields. Both correspond to a multistate Hamiltonian, 
where either the diagonal terms (chirping) or the off- 
diagonal terms (rf pulses) have explicit time dependence. 
In Ref . H the output coupling process was demonstrated 
experimentally for the sodium / = 1 situation, and the 
transition probabilities were in good agreement with the 
predictions of time-dependent three-state models. 

Both evaporative cooling and output coupling demon- 
strate the need to have analytic solutions for the time- 
dependent multistate models of the rf-induced dynamics. 
Although these models can be easily solved numerically, 
they are often used only as a part of a bigger theoret- 
ical description, in which the dependence of the solu- 
tions on the parameters such as the frequency and in- 
tensity of the rf field, or chirp parameters and shapes 
of pulse envelopes are required. Instead of looking into 
the known effects of rf fields, we can consider the ef- 
fects first, and then look how we need to tailor the rf 
field in order to achieve what we want — in this approach 
the analytic models show clearly their supremacy. Fur- 
thermore, as we show in this Letter, the description of 
the rf-induced multistate processes is closely connected 
to the two-state processes, which means that the wealth 
of knowledge on two-state models that has been accumu- 
lated in the past can be applied — and tested — with 
Bose-Einstein condensates. 

First we need to derive the Hamiltonian describing 
the rf-induced processes. The field B — Bq cos(w r f t)e r f 
couples to the atomic magnetic moment /i, i.e., H lnii — 

\jl ■ B. The matrix elements of this coupling between 
the magnetic states of the same hyperfine manifold are 
nonzero only if Am = ±1. Furthermore, using the angu- 
lar momentum algebra we see that the couplings between 
neighbouring states have the form [|]J|] 

H m , m +i = #m+i,m = yJJ- m)(f + m+ (1) 

where fl is the Rabi frequency quantifying the coupling. 
Here we have applied the rotating wave approximation 
and eliminated the field terms oscillating with frequency 
Lu r f . This leads to the curve crossing picture of the atomic 
potentials ||. 

In the regime of the linear Zeeman effect the energy 
difference between two neighbouring m states is Ez{R), 
which is independent of m but shares the i?-dependence 




FIG. 2. The curve crossing description of trap dynamics. 
The potentials are obtained by shifting the atomic states by 
multiples of rf photon energies. The region near the resonance 
point can be modelled with a bowtie crossing as shown in 
the enlargement of the grey area. Assuming that an atom 
traverses the crossing with some constant velocity v c , we can 
map the position-dependent crossing into a time-dependent 
one, which in the special case of two states corresponds to 
the widely used and well-known Landau-Zener model. In the 
optimal, i.e., adiabatic case the hot atoms follow the route 
marked by arrows out of the trap. 

of the trapping field. Here R is the distance from the 
trap center. Then 

H mm (R) = me(huj r t - \E Z (R)\) = mehA(R). (2) 

Thus A(i?) is the local detuning of the rf field. Since 
the trapping state can be either m = — f or m = /, 
depending on the particular atomic system, we need the 
factor e, which is +1 for m = —f trapping state, and — 1 
for m — f trapping state. 

In order to model the multistate dynamics we seek the 
solution of the Schrodinger equation for an iV-state sys- 
tem {Ti = 1): 

i-r-c = He, (3) 

dt w 

where c = (ci, C2, ■ ■ ■ , cm) t is the state vector containing 
the amplitudes for each spin state. For practical reasons 
we label the states with n = 1, 2, . . . N, instead of using 
the m labels (m = n — 1 — / = — /, — / + 1, ...,/; N = 
2/ + 1). The matrix elements of the model Hamiltonian 
are given by 
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H nn (t) = meA(t), 

H n , n+1 (t) = H n+l , n (t) = y/n(N - n)n(t), (4) 
= 0, (|n-fc| >2). 

Note that for a moving atom the -R-dependence in A 
can be mapped into time-dependence using a classical 
trajectory. 

We assume that initially the system is in the trap- 
ping state, which corresponds to either n — 1 or n = N , 
depending on e. However, due to the symmetry of the 
model we solve it for the case e — +1, i.e., start with the 
initial conditions 



ci,(-oo)=l, c n> i(-oo) = 0. 



(5) 



Then the case e = — 1 is obtained by reversing the state 
labelling and the sign of A. 

The model (|j) is a generalization of the Cook-Shore 
model, where A and O are time-independent |ic|j . We 
show that the solution for the iV-state model with the 
initial conditions (^) can be expressed using the solution 
(ox, a 2 ) of the two-state equations 
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Moreover, our derivation is considerably simpler and 
more straightforward than the one in Ref . [|l0| , where the 
Hamiltonian is diagonalised by means of rotation matri- 
ces using the underlying SU(2) symmetry of the model. 
We begin with N = 3. The Schrodinger equation is 

-A 

nV2 o nVz c. (7) 
o nV2 a 



d 
i—c 

dt 



We make the ansatz c\ — \\a\,C2 — A 2 aift 2 ,C3 = X^a^, 
substitute it in Eq. (Q), and obtain 

2iAidi = — AiAai + A2fiv2a2, 

iA 2 (dxa 2 + ax<22) = Q,^/2(\ia\ + A 3 a|), (8) 

2iA3a 2 = A 2 i7v / 2ai + AsAa2. 

By substituting cti and d 2 , found from the first and third 
equation, into the second we conclude that the latter will 
be satisfied identically if A2 = 2AxA 3 . Furthermore, it 
is readily seen that if we take Ax = A3 — 1,A2 = v2, 
the first and third equations for a\ and a 2 reduce ex- 
actly to Eqs. (||). Thus the solution to the three-state 
equations (R) is indeed expressed in terms of the solu- 
tion (ax, (12) of the two-state equations (||): c\ = af, C2 = 
\/2aia,2,C3 — a|- 

The result for TV = 3 encourages us to try in the case 
of general N the ansatz 



cx 

C2 



\ N-l 

\ n N ~ 2 



ft 2 



N-nn-1 



(9) 



and we choose Ax = Ajv = 1- We substitute this ansatz 
in Eq. (||) and from the first and the last equations we 
find the following equations for dx and d 2 



i(N - l)dx = -jAftx + A 2 V^ - lCla 2 , 
i(N - l)d 2 = Ajv-xV^V - max + jAa 2 . 



(10) 



By substituting these derivatives in the equation for c n , 
we conclude that the latter will be satisfied identically if 




Ajv-iA„, 



(11) 
(12) 



By changing n — > n + 1 in Eq. (11) and multiplying it 
with Eq. (ph we find that 



\ 2 \n-i=N-1 (13) 
By applying Eq. ( |l2] ) repeatedly n times, we obtain 



(N-l)\ 



1 



1)1 (N- 1) T 



(14) 



where we have accounted for Ax = 1. We now set n = N— 
2 in Eq. (|lj) , and taking Eq. (Q) into account we obtain 
A 2 = \n-i = VN — 1. Then Eq. (|l4| ) immediately gives 



An 



(N-l)\ 



(n - 1)\(N - n)V 



(15) 



Thus, we conclude that the solution to the TV-state equa- 
tions (||) is expressed in terms of the solution (ax, (X2) 
of the two-state equations (||) by Eqs. @ with A„ 
given by Eq. (plf). Furthermore, the A^-state initial 
conditions (||) require the twO-state initial conditions 
ftx(— 00) = l,a 2 (— 00) = 0. This implies that the final 
populations P n = |c n (+oo)| are expressed in terms of 
the two-state transition probability p = |a 2 (+oo)| = 
1 — |fti(+oo)| as 



P 2 = (N-l)(l-p) N ~ 2 p 

p - ( N - l V- (l - v )N-n n-l 

Pn=P N - 1 - 



(16) 



In estimating the spin-change probabilities for evapo- 
rative cooling we can use the Landau-Zener model pd| , 



A(t) = at, f2(t) = Oq = const, 



(17) 



where a is proportional to the change in A(i?) and to 
the speed of atoms, both evaluated at the trap edge (at 
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FIG. 3. The transition probabilities for the / = 2 multi- 
state Landau-Zener model. The two-state solution is given 
by p, and goes clearly faster to unity with increasing 0g/ a 
than the corresponding five-state probability P5. 



the rf resonance) . In the two-state model the transition 
probability is 



p = 1 — exp(— 2irti^/a) 



(18) 



Thus our model provides the exact result for the transi- 
tion probabilities in the generalized multistate Landau- 
Zener model, which can be used in estimating the ef- 
ficiency of the evaporative cooling M. As an example 
we show the / = 2 situation in Fig7]2|, where the final 
populations are plotted as a function of the adiabaticity 
parameter f^/a. As expected, in the multistate case one 
needs larger values of fl^/a to achieve population inver- 
sion, than in the two-state case. 

The Landau-Zener model can also describe the chirped 
output coupling. Then the atoms are assumed to be sta- 
tionary, so the time-dependence in A arises from the 
time-dependent change in cj r f (chirp), which is typically 
linear in time. In Ref. || the three-state version of the 
result ( |l6| ) was successfully used in describing the corre- 
sponding experiment in sodium / = 1 system. However, 
there the model was introduced only intuitively, and jus- 
tified merely by a comparison with numerical solutions 
to Eq. (J3J) . Here we have provided the proof that the 
/ = 1 solution is exact, and furthermore, derived the ex- 
act solution for any f . The special case of the three-state 
Landau-Zener model is also studied in Ref. fll2| . 

Instead of a chirped field one can use resonant pulsed 
output coupling, which was also demonstrated in Ref. || . 
For any resonant pulse we have A = and thus the two- 
state system follows the area theorem 



p = sin 2 (A/2), A = 2 



dttt{t), 



(19) 



duration). In the MIT experiment || the number of 
atoms left in the trap oscillated as a function of the area 
of a resonant square pulse, exactly as Eqs. 16 and jl9| ) 
predict. However, our model is not limited to resonant 
pulses only. For off-resonant pulses (A = const ^ 0) in 
two-state systems there are several known analytic solu- 
tions, which are reviewed e.g. in Refs. 

The purely time-dependent output coupler models de- 
scribed above are valid only if the time scales for the rf- 
induced interaction and the spatial dynamics of the con- 
densate are very different. In molecular systems one ex- 
pects interesting effects when the excitation process and 
internal dynamics of the molecule couple Q . It might be 
possible to realize some of the predicted molecular wave 
packet phenomena using condensates. 

In this Letter we have shown that any time-dependent 
multistate model describing the rf-induced coupling be- 
tween the different atomic spins states within the same 
hyperfine manifold can always be solved in terms of the 
solution of the corresponding two-state model. The fact 
that these models play a crucial role in time-dependent 
control of magnetically trapped ultracold atoms adds sig- 
nificantly to the importance of this result, which is also 
in general quite fascinating. 

This research has been supported by the Academy of 
Finland. K.-A. S. thanks Paul Julienne for enlightening 
discussions on evaporative cooling. 
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where A is defined as the pulse area, typically A tx £IqT 
(here fio is the pulse peak amplitude and T is the pulse 
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